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Abstract. We look at the time dependent fluctuations of the electrical charge in an 
open ID quantum system represented by a quantum dot experiencing random lateral 
motion. In essentially non-adiabatic settings we study both diffusive and ballistic 
(Levy) regimes of the barrier motion where the electric current as well as the net 
pumped electric charge experience random fluctuations over the static background. 
We show that in the large-time limit, t —> oo, the wavefunction is naturally separated 
into the Berry-phase component (resulting from the singular part of the wave amplitude 
in the co-moving frame) and the non-adiabatic correction (arising from fast oscillating, 
slow decaying tails of the same amplitude). Based on this separation we report two key 
results: Firstly, the disorder averaged wave function and current are asymptotically 
mainly determined by the same Berry phase contribution that applies in the case of 
adiabatic motion. Secondly, after a short transition period the pumped electric charge 
exhibits fluctuations that grow much faster than predicted by the adiabatic theory. We 
also derive the exact expressions for the average propagator (in the co-moving basis 
representation) for the diffusive and ballistic types of motion considered. 
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1. Introduction 

Since the pioneering work of Thonless [I] quantum pumping, i.e. transporting charge 
via a time-dependent potential or potential barrier (rather than applying voltage 
bias) has attracted great attention in the theory of mesoscopic devices, chiefly with 
the goal to mediate transport in Josephson junctions m ei m and quantum dots 
[3 ©00 El EB HUES]. The simplest theoretical framework for the quantum pumping 
phenomenon is provided by a ID quantum scattering at a time-dependent potential. The 
overwhelming majority of the existing results on time dependent quantum scattering 
were obtained assuming that either the time dependence of the potential barrier is due 
to periodic variations of its height and/or position, [9j EE} 'll] fl3lfITl To! EH} HJ or the 
changes are adiabatic so that at each time moment one can use a static scattering matrix 
©Eld to calculate the net transported charge. Some notable exceptions from these 
frameworks were considered in [18] where non-adiabatic pumping was considered in a 
discrete tight-binding model with harmonic on-site potential and [T9j where corrections 
to adiabaticity were obtained using a Kubo formula approach (see also Ref. [12] for the 
non-adiabatic spin pumping). Some of the numerical results reported in [11] also concern 
the non-adiabatical pumping but with AC variations of both height and position of the 
potential. Additionally, in Ref. [20] the impact of a randomly varying bond length in a 
quantum graph on the current was considered. 

In this paper we want to study quantum pumping in new settings where it results 
from non-stationary scattering at a narrow barrier (modeling quantum dot) experiencing 
a random walk in time, while preserving its shape. We consider a continuous time 
random walk (CTRW) where the scatterer can move either diffusively, via Brownian 
motion (BM), corresponding to delta correlated Gaussian velocity or with a constant 
velocity between the consecutive turning events with the “fat-tailed” distribution 
between these events (Levy walk) [211 , 22] , The latter case generally corresponds to 
superdiffusion. Both cases (especially the BM) provide bona fide examples of a non- 
adiabatic change of the parameters. Our goal here will be twofold: first, we will show 
that if one is interested in the transported charge then one can get a non-zero fluctuating 
value via lateral translation only (since periodicity is lost). This net charge should have 
zero average for unbiased CTRW and here we study how its variance grows with time. 

Our second goal is to show that in the non-adiabatic limit of CTRW an arbitrary 
initial state does evolve into the combination of a Berry contribution that follows the 
change of the scatterer position up to the same geometric phase factor [23] (see also 
0) and the non-adiabatic correction that leads to fast fluctuations of the probability 
current. The amazing result is however that when one considers the disorder-averaged 
values of the wave function and the probability current then asymptotically these seem to 
be solely determined by the Berry phase contribution. We call this result non-adiabatic 
Berry phase (BP) regime. In this asymptotic regime the average electric current is 
almost uniform and the growth of the fluctuations of the pumped charge is determined 
by the non-adiabatic fluctuations only. We show numerically that these fluctuations 
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grow faster than predicted by the commonly used adiabatic ” snow-plow” formula. As 
a by-product of our analysis we also obtain closed form expressions for the disorder- 
averaged quantum propagator (in the co-moving frame representation) for both models 
of the CTRW. Finally we note that although the role of the BP in the quantum pumping 
and quantum ratchets has been studied before [23, ESI ES EZ] (largely by considering 
some form of a quantum master equation) our paper presents the first results (to the 
best of our knowledge) on the strongly non-adiabatic stochastic quantum pumping in 
mesoscopic systems. 

2. Problem statement and the model 

We are considering a ID quantum dot placed between two leads at zero temperature 
at thermodynamic equilibrium. No external gate voltage is applied so both leads have 
the same value of chemical potential corresponding to the Fermi level: /i_ = fi + = Ep. 
Let us assume that the quantum dot experiences only lateral motion without any shape 
change of the potential. This can be induced e.g. by external AC irradiation of the 
scatterer. The situation can be modeled by the ID non-stationary Schrodinger equation 
with a given potential barrier V(x) experiencing a lateral shift in time: 

= + V ( X ~ 7 ^))^ W 

where 7 (t) is a time dependent position of the scatterer and we assume that 7(0) = 0 
always i.e. at t — 0 the system is fully characterized by the eigenstates of the static 
problem. We will assume here that the barrier is narrow compared to the typical wave¬ 
length of the electron wave (which is the Fermi wave-length, ‘lujkp). Without loss of 
generality we will consider a delta-shaped repulsive scatterer, although in Section [5] we 
also present some generic results valid for an arbitrary shape of the barrier. As for the 
position shift 7(t) we will assume that it experiences a CTRW of a particular kind. 
Specifically we assume that the barrier is moving with a piece-wise constant velocity 
and from time to time experiences a “jolt” which instantaneously changes the velocity 
to a new value after which the system continues to move with the new velocity until the 
next jolt, etc. We will assume that jolts represent events connected to a random process 
with independent stationary increments. Two different models are considered. In the 
first (Poissonian) model the time between the jolts (i.e. waiting time) is exponentially 
distributed and its average, r, plays the role of the velocity correlation time. As for the 
velocities we will assume that at each jolt the new velocity is chosen from a zero mean 
Gaussian distribution with the variance = 2D/r, with some positive constant D. 
The second case is the two state velocity model [21] (or Levy walk in the terminology of 
[22]) where the inter-event time has a heavy tailed distribution (i.e. decays as r a /t 1+a , 
0 < « < 1 as i 00) while the velocity at each event changes direction and can only 
take two values ±u. For the range of the exponents a given above all the momenta for 
the waiting time distribution starting from the first diverge and this process does not 
have a well defined characteristic scale m 
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The reason why these two particular models are chosen is that for the former model, 
in the limit r —> 0 (which will be studied here) one obtains the continuous white noise 
model for the velocity: 

t 

7(f) = f m dtW) = 0, = 2 D 8(t - t'). (2) 

o 

This corresponds to the unbiased diffusive BM with the familiar Einstein’s relation for 

the time growth of the root mean square displacement (RMS): yj 7 2 (t) = \/2Dt. On 

the other hand it is known that the Levy walk (in the specified range of the exponents 
a) is ballistic with the spread velocity V — (1 — a) 1 ^ 2 v. In all our numerical simulations 
we have chosen the so-called Levy-Smirnov distribution for the inter-event time, p(t), 
which corresponds to a = 1/2 [22j. Let us stress that we will assume that in both models 
the parameters of the CTRW are essentially non-adiabatic and the scatterer position 
experiences rapid random fluctuations. 

The solution if>(x,t) at an arbitrary moment of time can be presented as an 
expansion over Galilei shifted static eigenstates in the co-moving frame (see e.g. » 

/" iE(k) t , „ 

= I dkc(k,t) i>k{x - 7 (f)) e " , E(k) = —-. (3) 

The above equation defines a new time dependent shifted basis \k,^{t)) obtained from 
the eigenvectors of the static problem |/c, 0 ) via unitary transformation: 




U(t)\k, 0), 


U(t) = exp 



pj(t) ~ jH 0 t , 


(4) 


where Hq is the static part of the Hamiltonian in (JT]) and p is the ID momentum operator. 
In what follows we will mostly omit the y-dependence for brevity and use shorthand 
notation | k) for the co-moving basis. Clearly in the coordinate representation one has 
(x\k ) = *f>k(x ~ 7(t)) exp [—iE(k)t/h\. 

Of particular importance to the problem of the non-stationary quantum transport 
are the so-called scattering states Xk 0 [EE ESI- These are defined as the solutions of 
time dependent equation (JTJ) that at t — 0 coincide with the static eigenfunctions 
if>k 0 { x )- I n the co-moving frame representation (J3]) the scattering states correspond 
to the initial condition c(/c, 0) — 8(k — kf). Another interpretation of a scattering 
state Xk 0 is that when projected on the moving basis | k) it provides the expression for 
quantum propagator in this representation (see Section [5]). Therefore these states can 
also be used to study random scattering of narrow wave-packets in the k- space. When 
the potential experiences purely periodic motion starting from t = — oo the scattering 
states can be obtained in the closed form via the Floquet formalism PHOj. When the 
time-dependent potential is small one can also obtain these states perturbatively |28]. 
In general however these solutions must be obtained either numerically or by making 
some additional assumptions e.g. considering the asymptotic limit t —> +oo which is the 
method we adopt in Section [4j An important limiting case is that of the adiabatically 
slow motion of the barrier. When the change of position is adiabatic then at time t the 











5 


scattering state is given by the corresponding moving eigenstate \k 0 ) np to a BP factor 

mm 


|x*o(0) = e^‘oW|fc„). 


( 5 ) 


So our first goal in this paper will be to see what happens to the amplitudes c(k, t ) of 
the scattering state Xk 0 { x ,t) when the changes of the potential are fast and random as 
assumed in our model. As we shall see, in this case the complex amplitude c(k,t ) is 
characterized by sharp delta-peaks at k — ±k 0 while at large times the wave function 
in the coordinate representation, Xk 0 ( x , t) naturally separates into the local Berry phase 
contribution ([5j) and nonadiabatic corrections arising due to the slow (power law) decay 
of the amplitude c(k,t ) at \k\ — y oo. For the Levy model the situation is more 
complicated since the system spends relatively large times moving with constant velocity 
±n 0 which lead to additional Doppler shifted peaks in \c(k,t)\. 

We will also be interested in the properties of random quantum pumping in such 
a system and more specifically in the pumped charge. Traditionally when considering 
adiabatically slow periodic changes of the parameters of the scattering channel (i.e. 
scattering matrix) the pumped charge is calculated via the well known formula of 
Biittiker, Thomas and Pretre (BTP) [5]. For a single channel system where only lateral 
shifts of the scatterer are considered (i.e. our situation) this formula can be considerably 
simplified giving for the infinitesimally small pumped charge mm 

AQ = -\r kF \ 2 Ay, (6) 

7T 

where r kF is the static value of the reflection coefficient at Fermi level and Ay is the 
infinitesimally small change of the scatterer position. Eq. () 6 |) reflects the so-called “snow¬ 
plow” contribution to the charge [7j. 

Since our paper concerns with the non-adiabatic and non-periodic quantum 
pumping one will need to simulate the electric current passing through the leads 
directly, without simplifying assumptions. Using a standard procedure in which different 
scattering states serve as a basis for expanding the time and space dependent held 
operator T (x, t) in terms of creation and annihilation operators of the electron in the 
right and left lead and performing quantum mechanical averaging at zero temperature 
one arrives at the following expression for the charge pumped through the left or 
right “+” lead [28], [30] determined at the positions ±L/ 2: 

t kp 

Q±(t) = Jl±(t')dt ', I±{t) = -e J J ko (±L/2,t)dk 0 (7) 

o -k F 


where e > 0 is the elementary charge and 


is the probability current (flux) associated with the particular scattering state Xk 0 - 
The suitable choice of the boundaries ±L /2 will be to place them outside the RMS 


Jk o i x i t') 


h 


2m i 


X* ko ( x ,t) 


dXk 0 ( x ,t) 

dx 


Xk 0 ( x ,t ) 


dxl 0 { x ,t) 

dx 
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displacement of the random potential. The physical meaning of Eq. (17|) is transparent: 
the electric current passing through right (left) lead at time t is proportional to the 
total probability flux of all right ko > 0 and left ko < 0 moving electrons under the ID 
“Fermi surface” [— k F , k F \. For the static case 7 = 0 the scattering states are just phase 
rotated eigenstates \Xk 0 ) — so that the left and right moving currents are 

proportional to the total transmission for the right and left movers (the celebrated 
Landauer-Biittiker formula [30], [31, [32]). In the absence of gate voltage, for any 
symmetric potential, the time reversal symmetry ensures that these two contributions 
exactly cancel each other and the total current in each lead is exactly zero (see e.g. [30]). 
Introducing time dependent motion of the quantum dot generally breaks time reversal 
symmetry, however for adiabatically slow periodic pumping the net charge accumulated 
over one period of oscillations is still zero m- In our case of random non-adiabatic 
pumping the total charge Q(t ) represents a random process which is expected to have 
zero mean and time dependent RMS <7<g(t) which is the quantity we are interested in 
here. 

Since the scattering states form the basis of the theoretical treatment of the non- 
adiabatic quantum pumping in the next few sections we will study in detail the statistical 
and asymptotic properties of these states at arbitrary wavelength A = 2n/k 0 of the 
incoming wave. The results of the following three chapters are quite generic and can 
be applied beyond the context of quantum transport in such areas as non-stationary 
scattering of wavepackets HZ] or quantum graphs with fluctuating bond lengths |20j . 


3. The dynamics of scattering states 


Since each scattering state is a solution of the time dependent Schrodinger equation (JTJ) 
it is natural to expand it in the co-moving basis (|3]). Let us now derive the continuous 
time evolution equations for the amplitudes in the co-moving frame c(k,t). Inserting 
expansion (]3]) into the Schrodinger equation (JT|) and projecting it onto the corresponding 
moving eigenstate one obtains the following exact equation [13] : 

c(k : t)= A f(t ) J A(k,k'',t) c{k' ,t) dk ', c(k, 0 ) = 8{k — ko). (9) 

Here the time-dependent kernel A(k, k'\ t ) is given by 

A(k, k'-1) = e* dRD-WR f ^*( x ) ^ ^,( x ) dx = (k, 0 | e i6 ° t/h (ip/h)e- i6ot / A \k', 0 ). ( 10 ) 

J ox 

Note that this kernel depends only on the static Hamiltonian. Also it is interesting that 
Eq.flnD with the kernel (HUD looks like the evolution of the wave-function in the interaction 
(Dirac) representation (see e.g. [33j) with the role of perturbative interaction played by 
the operator 7 (t)ip/h 2 . 

We are using the model of a delta-barrier V(x) = ( h 2 /m ) fl<5(a;) with a continuous 
twice degenerate spectrum. Since the barrier is symmetric the static eigenfunctions 
corresponding to the left traveling waves k < 0 are obtained from the right-traveling 
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ones by the coordinate flip: ^-fc( x ) = The scattering is completely described 

by the unitary S— matrix: 


S = 


t' 


Q 


Tk = 


ik — 0 ’ 


4 — 


i k 


i k — fi ’ 


4 = 4 , r' k = r k { 11 ) 


and the left/right traveling normalized eigenfunctions are given by 

M x ) = ( 2 vr )- 1/2 {t w e ikx 6[±x\ + ( e ikx + r ]k , e~ ikx ) 0[ T x]} , 

where the signs of the arguments of the Heaviside function 6[x\ are chosen according to 
the sign of k (reflecting the mirror symmetry). From the above one readily obtains an 
expression for the kernel A(k, k!\ 0 ): 

U\ _ " " ' k‘k'/ I 

'V o-r/’k/_M I 


ik \tk\ S(k > 


k \ t t t k>- 1 - r k r k>) , k '( r t+ r k') 


A(k, k 0 ) = < 


ik \t k \ 2 5{k! — k) + 


2'K(k'—k) 

k'(tkt* kl -l-r k r* k ,) 


2tt( k'-\-k) 5 


^1 u 1 27v(k , —k) 2 tt (k'-\-k) 

ih (+, r* - \ s(v j-tu k '( r *k t V- r t> t V M 

ik 
2 


— T (4^ - tfcr fc ) + fc) + 

-f (4 r, fc-44)h(4 + /c) + 


27r(/c , +Ai) 

k'(r k ,t k -r k ty) 


+ 

+ 


27r (k'—k) 5 
k'jtk-ty) 


k > 0 , k 1 > 0 
/c < 0 , k! < 0 
k > 0 , k' < 0 
k < 0 , > 0 


.( 12 ) 


27r(/c , +Ai) 27 t (k'—k) ’ 

Generally for any symmetric potential it follows that the kernel is asymmetric with 
respect to inversion: A(—k,—k']t) = —A(k,k']t). This means that if we know the 
solution of Eq.Q, say c(k, t ), then c(—k, t ) is the solution of the time-reversed equation 
Eq. fj9l) with 7 (t) —> —j(t). Also since the Hamiltonian H 0 generally does not commute 
with the momentum p the matrix A(t) is always non-diagonal. 

Another technical note is that if y(t) is the Brownian motion, (l2]h then Eq. (191) 
represents a stochastic integro-differential equation with multiplicative white noise. This 
means that one can understand it either in Ito or in Stratonovich sense [33]. Here we 
will opt for the latter, which actually implies that the two consecutive jolts separated by 
infinitesimally small average time moments are correlated - which seemingly contradicts 
the independence assumption postulated above in our model. We argue however that 
this discrepancy (expressed in the additional spurious drift [31]) is not crucial for the 
observed physical effects and is simply an artifact of the chosen model. Indeed the latter 
can be always amended to include correlated velocity increment to yield the Stratonovich 
picture in the limit r —> 0 . 


4. The asymptotic Berry phase 

In this section we shall mostly concentrate on the asymptotic behaviour of the solution 
c(k, t ) of (jS]) for a given realization of the CTRW q(t). The ensemble averaged scattering 
states are given in the next section. Many of the results of this section are valid for 
arbitrary continuous translations of the barrier and may have applications beyond the 
scope of quantum transport. 

For an arbitrary barrier shape it is not possible to obtain the solution of flU]) in a 
closed form. However an important simplification comes from considering the structure 
of the time-dependent kernel A(k, kt) in Eqs. (1101120 . Indeed one can see that the first 
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two terms in A(k, k 0 ) always contain singularities at k — ±k' (the delta-like singularity 
and a simple pole) so that the main contribution of these terms in the r.li.s. of Eq.Q 
comes from the vicinity of these poles. In particular if one starts with the pure (right- 
moving) eigenstate |/c 0 )> than instantly the opposite, mirror-reflected eigenstate | — k 0 ) 
is excited. These states continue to move in space, dragged by the barrier, so that at 
any given time moment the amplitude in the co-moving frame c(k, t) will contain the 
two peaks of variable height at k — ±k 0 . But in the vicinity of these peaks one can 
neglect the oscillating term exp [i(E(k) — E(k'))t/H] in the definition of the kernel A 
(HO i.e. assume that A(t ) ps A(0 ). This approximation fails to describe the tails of 
the amplitude c(/c, t) and, as we shall see below in the limit of large time these tails are 
responsible for the break-up of the adiabatic BP ansatz. 

For the static kernel A(0) the solution of Eq.(j9]) is expressed via the usual (not 
time ordered) exponential propagator exp[ 7 (f)A( 0 )]. As mentioned before the operator 
A( 0 ) = ip/h is not diagonal in the |fc, 0 )-basis since it generally does not commute 
with the Hamiltonian. This difficulty can of course be overcome by changing into the 
continuous eigensystem |A) of the momentum operator which is just a set of normalized 
exponentials exp [i Ax] . In this basis we can easily calculate the propagator analytically 
and then return into the original k representation to get the solution of Eq. (l9l) . 

It is straightforward to show that if initially the system is in the eigenstate \k 0 ), 
then the amplitude at the time moment t is given by: 

c(M) = (k\t) = I d\e iXl{t) (k\\) (A|fc 0 >, (13) 

where real line singularities (if any) should be treated in the sense of Cauchy principal 
value. 

Next using the static wave-functions of our delta-barrier in the coordinate domain, 
tjjkix) one obtains: 

/ I X \fv T** 1 "I 

dx Vk{ x ) e * Xx = ^ A 2 Jp) + 2^ + ^fc|) 5 ( A - k ) + 2 r \k\ 5 ( X + k ) 

(ai k ) = m *. 


Using the expression above and performing Fourier transform with respect to A fEq. (fT3l) l 
we finally arrive at the following answer: 

c(k, t) = f(k 0 , t) 8(k — fc 0 ) + g(~ k 0 i t ) 5(k + k 0 ) 


2n (k 2 —k‘o) 

i k O r ko 


2 Sign[/c] r* fe| r feo [k sin(/c 0 | 7 (t)|) - k 0 sin(fc| 7 (t)|)] 


(! + t* k \ + r |fc|) cos(fc7(t)) + i (1 + - rjy sin(fc7(t)) 

-i \k\r\ k \ [(1 + t ko + r ko ) cos(k 0 j(t)) +i(l + t ko - r ko ) sin(fco 7 (t))]} , 
where in order to obtain the expression above we have made use of the identity: 


(14) 


7 A 7 


dX 


l-oo (A 2 -F)(A2-^) 

7 r /csin(/c 0 | 7 |) — A: 0 sin( A: I 7 1) 


k kn 


U 2 _ U2 

ft A/Q 


7 r 

+ ^2 _ k o) + 5{k + k 0 )) cos k 0 7 (f), 


p.v. 
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where the regular part of the expression above can be easily obtained via residues and 
the delta-part (corresponding to coalescing poles) can be calculated by integrating the 
l.h.s. over two infinitesimal semi-infinite stripes \k ± ko\ < e, interchanging the order of 
k and A integration and letting e —> 0. 

The amplitudes f(k 0 ,t ) and g(—k 0 ,t ) have simple closed-form expressions (for a 
delta potential): 

f{k 0 ,t) = cos(fc 0 7 (f)) + i + sin(/c 0 7 (t)), g(-k 0 ,t ) = - ^^2 sin ( fc o7(*))■ (15) 

and correspond to unnormalized Weyl differentials [22J, [SB] in the moving basis 
representation of the initial state and its mirror counterpart. 

One can see that the result (TT4h . (1T5T) is quite far from the adiabatic BP formula (J5l) 
as the states different from |fc 0 ) are excited and the amplitude of the state itself, f(ko, t ), 
(which has a natural meaning of survival amplitude) is not a pure phase factor (unless 
the scatterer is very weak, 0 <C ko). But as we shall see below, when one considers the 
wave function in the coordinate domain the result ([5]) is largely recovered. 

At this stage it is pertinent to see how well the results above describe the genuine 
behaviour of the amplitude c(k,t). To this end we now turn to numerical simulations 
performed within the framework of the delta-scattering potential model (fTTlh(fT2lh We 
have solved Eq.([9]) assuming that fc-space is discretized in the range [—3/c 0 , 3/c 0 ] using 
4096 discrete points with the grid size, A k, imposing an effective box quantization 
with the box size L' = 2ir/Ak. The spatial coordinates were measured in units of 
l/k 0 while time was measured in units of the wave-period T = 2nh/E(k 0 ). The 
amplitudes themselves were also rescaled: c(/c, t ) —)■ c(/c, t ) A k so that the initial 
condition c(/c, 0) = 5{k — k 0 ) now corresponds to a single discrete peak ko . In all our 
simulations we have chosen the values of the parameters that correspond to the non- 
perturbative regime and the scattering was “strong”: hi = 0.5fc 0 - As for the specifics of 
the CTRW we used (in the selected units) the values D = 0.5 for the BM while for the 
Levy walk the constant velocity vq was chosen in such a way that for the corresponding 
wave-vector k = mv^/Ti one has k = 0.3fco- For the Levy walk the velocity changed 
sign after random time moments sampled from the Levy-Smirnov probability density 
function (PDF) P(t ) = r 1 / 2 /(2 A /¥t 3 / 2 ) exp(—r/4t) with r = 4x 10 -3 (which is the same 
as the smallest time grid in the diffusive case). 

The results for the amplitudes c(/c, t ) and f(ko, t ) are given in Fig]T]for two separate 
realizations of the Brownian motion and Levy walk. The first result of the comparison 
is that Eq. (lT4l) indeed provides a very good approximation in the vicinity of the peaks 
k = ±ko. To demonstrate it in FigOJb) we have enlarged the area aroimd the main 
peak k = k 0 . Both real and imaginary part of /(/co,t) are in excellent agreement with 
the predictions of (TT5|) as seen from Fig{l](c),(d). 

However outside the peaks the difference c(k, t) between the genuine amplitude 
c(k,t) and the approximation co(k,t ) provided by Eq. (lT4l) is obvious. One can see 
that c(k, t ) is a regular, fast oscillating function with slow decaying power-law tails. 
Another feature concerns the difference between Brownian and Levy motion. First 
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■ 0.1 

, 0.01 


\ 

-3 -2-10 1 


-BM simulations 

-BM theory 

Levy simulations 
- Levy theory 


k/k 0 




0.92 0.96 1.00 1.04 1.08 
k/k 0 



Figure 1. The comparison of the modulus of the numerical amplitude \c(k,t)\Ak 
after t = 8 wave periods, (a),(b) and the survival amplitude f{ko,t ), (c),(d), with 
Eqs.(fni)-(fl5l). 


of all the Levy case is characterized by additional Doppler shifted peaks located at 
k = ±|fco ± 2A;|. These peaks characterize the motion with constant velocity ±t>o and 
are discussed in more detail in section [5] (see also Ref. EH])- The most significant of 
these peaks are denoted by arrows in FigOJa). Note that since for the BM case the 
average instantaneous velocity (or k ) is zero these peaks are absent. Also since in the 
Levy case the system spends a relatively large amount of time locked into a constant 
velocity motion a large part of the corresponding graphs [0(c), (d) displays just sinusoidal 
oscillations as follows from (TT5]l . 



Figure 2. The integration contour for evaluating the wave-function Q via m- 

Next let us turn to the time evolution of the scattering state in the coordinate 
representation as it is required to define the probability current (J8]). To do so one must 
plug the amplitude (fIT|l into expression ([2]) (understood in the Cauchy principal value 
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sense). The delta-terms obviously do not present any difficulties while in order to deal 
with the real axis poles k = ±/c 0 we choose the loop in such a way that the exponential 
exp [—iE(k)t/h] decays everywhere apart from the real axis and all the poles of the 
transmission and reflection coefficients are avoided. We have opted for the contour 
shown in Fig 12 In the limit of large t the integral (j3J) over the bisector 9ft [k] = — A [4] 
is of the saddle point type with the main contribution coming from the origin and it is 
straightforward to show that it decays in time as t ~ 3 / 2 (see Appendix A). On the other 
hand the residue contribution in (1T41) combined with the delta terms provides a simple 
asymptotic result for the scattering state: 


Xk 0 {x, t) = - (1 + 9ft(4 0 - r k o - r ko t* ko )) e lkol{t) 'i/j ko (x - 7 (t)) e 


iB(k n ) 
h 


— e *feo7 (t) 


/ / \ iI2(fco )f 

^koix -7(t))e t y 00 . 


(16) 


But this looks remarkably like the Berry phase result ()5]) in the coordinate representation 
obtained here by taking into account only the singular part of the wave amplitude c(k , t). 
The effective Berry phase here is <p ko (t) = /c 0 7 (t). The same result one gets by applying 
the formula for the BP for a continuous eigenstate [29] which in our notations reads 
0fco — il{t) f A(k, k 0 ] 0)dk 0. In the same approximation the probability flux built on 
the wavefunction (|T 6 |) is just the static probability current of the initial eigenstate 4> : 

J ko (x,t) = J s t(k 0 ) = hk 0 \t ko \ 2 /( 2 nm), t ->• 00 . (17) 

We will call the results (Il6|) . (fl7|) the Berry phase contribution. Recall that these were 
obtained in the limit t —> 00 by neglecting the oscillating tails of the amplitude c(k,t). 
To see the nature of these fast oscillating tails let us return to Eq. (l9]l and present the 
solution as a sum c(k,t) = c 0 (k,t) + c(k,t) where c 0 (k,t) is given by Eq. flTT]) and is the 
solution of (jH]) if one neglects the time dependence of the kernel A. Then for the 11011- 
adiabatic correction c(k,t ) that we assume due to numerical evidence to be a smooth 
differentiable function one obtains an inhomogeneous equation: 

c(k,t) — j(t) j A(k, k'\t) c{k!,t) dk! = 7 (f) J [A(k, k'] t) — A(k, A/; 0)] c 0 (k', t) dk'. (18) 

We can evaluate the right hand side exactly using the expressions for A(k, k',t ) in Eq. 
(fT2]i and for co(k,t) in Eq. (fT 4 lh As we assume furthermore that c(k,t ) possesses no 
singularities, we perform on the left hand side the remaining k'- integration for large t 
by closing the integration contour in the complex plane using the same integration loop 
as in Fig]2 We then obtain the following equations with k > 0 

c(k, t ) - 7 (t)ik (| 4 | 2 - \r k | 2 ) c(k, t ) + 7 (t)ik ( t k r* k - r k t * k ) c(-k, t ) (19) 

= 7 (t) ( e i(m-E(ko))t/H F ( fc> jfco> e iky(t) 7 e <*07(*)) + G(k, k 0 , , e iko7(t) )) 

and 

c(-k, t ) + 7 (t)ik (|4| 2 - \r k | 2 ) c(k, t ) + 7 (t)ik ( t k r* k - r k t* k ) c(k, t ) ( 20 ) 

f Strictly speaking applying the corresponding formula from [29] to the kernel (fl2ll gives rise also to an 
imaginary part of the phase. The reson for this is that in their derivation the authors of [29] assumed 
that the kernel A(k, k'; 0) is always diagonal which is clearly not the case for (1121) . 
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= j(t) (e i{E ^)- E ^))t/np ^ k? ^ e ifc7(t) > e ^o7W) + ^(fc, fc Q> e <M*) ; e <*o7(*)) j 

with the functions F(fc, jfe 0 , e ik ^\ e ik ° 7 ^), G(fc, fc 0 , e ifc7(t) , e ifco7(t) ) and F(fc, fc 0 , e**^, e 4 * 07 ^), 
G(fc, fco, e lkl ^\ e* fco7 ^) obtained by integrating the right hand of Eq. f[T 8 ft . We assumed 
that the function c(k , £) possesses no poles arbitrarily close to the real axis. Note in 
this context that the expressions on the right hand side of Eq. (1T8li and Eqs. (1191201) 
also have no poles in that region as the poles in co(k',t) are canceled by the factor 
e %{E(k)—E{k'))t/h _ i j n (JJgj) The system of differential equations shown above can 
be diagonalized and the resulting differential equations can be solved in the case that 
averages of functions of 7 (£) are known. The solution can be represented as time integral 
with the integrand also containing the inhomogeneities on the right hand side of Eqs. 
(U9l |20[) . Due to their time dependence oc e l ( E ( k )- E ( k o))t/n a p , 0 g(/c,f) will oscillate in 
that way in case there are no significant oscillations from other source, i.e. the noise. 



BM, simulations 



k 0 x 


-1 0 1 
k 0 x 


Figure 3. (a) The simulated scattered state |yfc 0 (x. t)| compared to the BP 

contribution, (b) the normalized probability current at t = 8. The deviation of the 
latter from the static value, 03 are due to non-adiabatic effects. 

All this is illustrated in Figj3] where the BP contribution to the scattering state (fT 6 l) 
and its probability flux (TTTT) are compared with the results of numerical simulations with 
the same parameters used to obtain FigJTl Only part of the BM solution is shown. In 
Ha) one can immediately discern the final position of the barrier (indicated by arrows) as 
the shape of the solution is different at the preceding and trailing ends. One notices that 
the non-adiabatic contribution, c, is important leading in particular to the non-uniform 
distribution of the probability current in[3](b) (note however that the fluctuations appear 
to be smaller for the Levy case). 

The latter result should not come as a surprise, after all one should not expect that 
the Berry phase result obtained adiabatically should hold exactly for a rapidly changing 
position of the barrier. However in the next section we shall see that surprisingly when 
one considers ensemble averaged quantities that are linear in the amplitude c(k,t), in 
the large time limit Eqs. (fl4|) and (H71) can be used to determine the average quantities 
while the non-adiabatic correction seems to be always averaged out. 
































13 


5. The disorder-averaged quantum propagators 

As mentioned earlier one important property of the scattering states apart from their 
being used as the basis for calculating the electrical current (J7J) is that expressed in 
the co-moving frame representation \k,^f(t)) (|3]),(j4|) they serve as quantum propagators. 
By these we means the solutions, K(k,k'-,t), of Eq. (f9]) given the initial condition 
K(k,k']0) = 5{k — k'). In the previous section we have already seen that for each 
realization of the CTRW each of these solution separates into the BP contribution (1141) 
and the non-adiabatic part. In this section we show how one can calculate the exact 
average of the propagator without any a-priori assumptions on the structure of the 
solution. The methods used here differ slightly for the BM and Levy case but before 
we consider each case in detail a few remarks are in order. Firstly according to the 
time-reversal symmetry mentioned earlier it is only sufficient to consider the averaged 
propagator K(k,k']t ) for positive k' as for each realization of 7 {t) the solution for 
the negative k! is obtained from its positive k' counterpart by reversing the direction 
of motion: k! — > fc, 7 — > —7. Since we here only consider a symmetric (unbiased) 
CTRW the average of any function of 7(t) is insensitive to the sign so one gets 
K(k, —k'\t) = K{k , k'\ t). 

Another important feature shared by all such types of CTRW is that if one is 
only interested in the average survival amplitude (or Weyl differential) of the initial 
eigenstate |fc 0 ) then as we have seen from Fig{l](c,d), Eqs. dlSj) are in excellent agreement 
with the numerics for each realization. Therefore one can average these formulae 
directly and arrive at the following amazing statement: Regardless of the nature of the 
symmetric CTRW of a narrow potential barrier the average survival amplitude f(k 0} t) 
of an eigenstate \k 0 ) is just the characteristic function &(k 0 ,t) of the current position of 
the scatterer 7(f). This result is more amazing inasmuch as it coincides with the disorder 
averaged BP formula ([5]) but obtained here in the non-adiabatic regime. In other words: 
if the system experiences symmetric CTRW the celebrated BP result is recovered on 
average. The characteristic function of the BM is readily obtained since 7 (t) in this 
case is a zero mean Gaussian variate with the variance 2 Dt. The case of Levy walk is 
slightly more complicated but one can derive ( |22] , chapter 8) an asymptote of the time 
Laplace transform <L(s, ko) for the case of the Levy-Smirnov distribution (cc = 1/2) for 
small s and v 0 k 0 r <C 1: 4>(s, k 0 ) ps (s 2 + v^kf)" 1 ^ 2 . The inverse Laplace transform of 
the above is known and provides the asymptotic form of the characteristic function as 
t —> 00. Summarizing the results one has for the average survival amplitude: 

ko+8k 

f(k 0 ,t)= lim / c(k,t) dk = 

Sk^-0 J 

ko~8k 

Note that the average amplitude is purely real owing to the symmetry of the CTRW. 

To see how well these results are confirmed by the numerics we present a result of 
averaging over 1000 realizations of the CTRW with the rest of the parameters identical 
to those of the previous section. The results are given in FigJdJ One can see an excellent 


g-Dfcg t 


BM 


J 0 (k 0 v 0 t), t^> t Levy (-Smirnov) 


( 21 ) 
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agreement between theory and numerics and observe the qualitative difference of the 
amplitude behaviour: the BM amplitude decays exponentially in time while the Levy 
walk corresponds to slowly (oc 1/t 1 / 2 ) decaying oscillations. The imaginary part of the 
average amplitude (not shown) was observed to be zero up to the third digital place. 

So how can one explain that the behaviour at ko is to a good approximation 
unaffected by the transitions to other /c’s? One way to understand it intuitively goes 
as follows: Starting out at t = 0 we have c(/c, 0) = S(k — k 0 ), i.e. the only excited 
state is |/c 0 ). From Eq. (f9j) we see that the initial evolution of this state’s amplitude 
is given by the diagonal (singular) part of the kernel, A[k,k',t ) which after averaging 
around ko, Eq. ()2Tl) . gives exactly the behaviour prescribed by the Berry formula (J5]). In 
the non-adiabatic case the other states are excited very quickly but since the CTRW is 
symmetric (i.e. unbiased) the average total amplitude transfer from the original ko to 
other /c’s and back is almost zero and the main effect again comes from the diagonal 
“self-action” of the amplitude c(k 0 ,t). For the BM case this can be rigorously justified 
which is done in Section f5T2l 



Time, wave periods 

Figure 4. The dynamics of the average survival amplitude compared to theoretical 
prediction. m 

Besides the averaged survival amplitude also the disorder averaged probability 
current is of interest. The results for the current are given in Fig. 0 The averaging was 
performed under the same conditions as in the previous sections: 1000 realizations of 
CTRW were used with the same parameters for each realization as those in the preceding 
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figures. For both BM and Levy CTRW the average area explored by the potential is 
shaded in grey. We can see that for the case of BM the averaged probability current is 
very close to its uniform static value which agrees well with the BP result (ITT]) . On the 
other hand for the Levy case the non-adiabatic contributions survive even after averaging 
leading to small oscillations of the average current. This is to be expected since in the 
Levy case the additional Doppler peaks of the amplitude c(/c; t) (see Fig. [0(a)) are not 
taken into account by Eqs. o and (1161) and their nonadiabatic contribution leads to 
the observed spatial oscillations of the average current. 


(a) (b) 



k 0 x 


Figure 5. Average probability current normalized to its static value J st (fco) after t = 8 
wave periods, (a)BM case, (b)Levy walk. 


Let us now turn to the averaged propagator and evaluate not only the BP 
contribution but also the average non-adiabatic diffusion in the fc-space. We start with 
the case of the Levy walk. 


5.1. Levy walk 


Since in the Levy model the scatterer moves freely between the two jolts for relatively 
long times it is natural to study the propagator for the barrier moving with constant 
velocity v, K v . In the coordinate representation one can of course use the Galilei 
transform of a static propagator K st (x, x'\ t ) of Eq. (JT]) for 7 = 0. The result is 

K v (x,x'-t) = e~ iECk)t/h K st (x - vt,x']t) , (22) 


where k = mv/h is the wave-vector corresponding to the velocity v. On the other hand 
the explicit expression for the static propagator is known for some important classes of 
barriers, including delta- and hyperbolic secant ones [35] [36]. For the time evolution of 
the initial state |/cq,0) in the coordinate representation we then obtain from (]22l) 


Xk 0 (x,t)= I dx'K v (x, x'; t)ijj k0 (x') 

dx'dke-^^^^-^^x - vt)^* k {x')^ kQ (x'), (23) 


' —OO 
/•OO 


















16 


where we have applied the eigenfunction expansion in the second step. We can perform 
the x '-integration exactly using the explicit expressions for the wave-functions of a delta- 
scatterer and afterwards evaluate the fc-integrals in the limit of large t using the contour 
in Fig. [2] to obtain (fc 0 , k > 0 is assumed): 


Xk 0 {x,t) 




/ _„ik 0 x ( ,—iE(k 0 )t/h 

iE(ko)t/h | ty, ~ V 2 


X > vt 

iE(k 0 -2k)t/h x < vt 


t —> OO. 


(24) 


The wavefunction corresponding to the barrier moving in the opposite direction is 
obtained via flipping the sign of the corresponding wave-vector k. Focusing on the 
transmitted wave or current at x > 0 we observe that the transmission coefficient t ko 
obtained for a static potential is now replaced by t ko _ k and t ko+k which inevitably leads to 
the appearance of electric current. It is interesting to observe how in the adiabatic limit 
of small velocity k <C /c_f one recovers the “snow-plow” formula ()6]) (see Appendix B ). 

When one considers a Levy walk however it is more advantageous to work with the 
expressions for the constant velocity propagator in the co-moving basis (J3]), (J4J) : 


Appendix B 


K v (k,k'-t) 


- E(k) t 
= e h 


J dk''(k\e <ii \k")ei lEih) - E{h " )n (k"\e- ai \k'). 


(25) 


The easiest way to derive the expression above is to notice that an arbitrary solution 
of Eq. (JT]) with a uniformly moving barrier can be presented as a linear combination 
of the moving frame solutions 0fc(x — vt) exp (ikx) and use the initial condition 
-0(^,0) = to find the coefficients. 

The propagator (l25]i has some features interesting for the future analysis. First 
of all matrix elements (k\ exp(±i kx)\k') can be calculated explicitly for a delta-barrier 
and they have delta peaks (as well as simple poles) at the Doppler-shifted wave vectors: 
\k ± k'\ = k. Thus the amplitude c(k,t ) initially corresponding to an eigenstate |fc 0 ) 
(/co > 0) will have six sharp peaks. Apart from the two main peaks at k — ±ko 
there are four additional ones located symmetrically at the Doppler shifted points 
k = ±(fco + 2/c) and k = ±|2 k — ko\. The two main peaks k = ±ko (always 
dominating in terms of magnitude) experience periodic time oscillations with the 
frequency uj± = (1 /h)(E(k) + E(ko) — E(ko ± k )) = ±k 0 v which coincides with the BP 
result ([5]). This is also consistent with the long time asymptotic result in the coordinate 
domain (l24lh The side peaks however cannot be obtained from a simple adiabatical 
approach of e.g. [7j and the self-consistent adiabatic treatment must include the newly 
created Doppler shifted states as the additional source of pumping for the main mode 
k 0 = kF (see EH]). The beating frequencies of the Doppler-shifted peaks are given 
by the combinations c o = (1 /h)(E(k) + E(k Q ± 2k) — E(k 0 ± k)). Note however that 
unlike Ref. m here these peaks are present in completely non-adiabatic settings. All 
this is illustrated in a specific example in Figj6]in which we considered a right moving 
scatterer with the same constant velocity as in the Levy walk example in the section 
(J4]) corresponding to k = 0.3/co- Most of the peaks are clearly visible and indicated by 
arrows. Note that not all six peaks can be observed at any given moment, e.g. the 
leftmost peak at k = —1.6/cq hr Fig. [6] is absent. 
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10' 1 : 



10 ‘ 3 - 


10' 4 

-2 




k/k 0 

Figure 6. The amplitude distribution for a uniformly moving scatterer after t = 8 
periods of the moving wave 



This picture helps to shed a new light on discrepancies between the /c-space diffusion 
for BM and Levy model presented earlier in FigOJa). Indeed for the Brownian motion 
the velocity of a particle at each time moment is sampled from a symmetric Gaussian 
distribution which has zero mean and a very large variance (inversely proportional to 
the correlation time, r). Therefore the effective Doppler wave vector k experiences large 
fluctuations around zero at each time moment and therefore no coherent picture of the 
peaks emerges. The Levy walk on the other hand spends a long time locked in the 
motion with fixed positive or negative values of the Doppler shift ±k which explains 
why at least some of the peaks from Fig. |6] survive in the stochastic motion in FigOJa). 

Next let us turn to the averaging of the multiple collision propagator, originating 
from a large number of velocity jumps. For our Levy walk it is the propagator K v 
(either in coordinate representation (122 ] ) or in the moving k- frame, (l3ll25|) l explicitly 
depending on the time interval and velocity, that serves as the building block for the 
overall propagator after n discrete random jolts: 

n 

K n = K v (Ati , v{), K 0 = I- 

i =0 

Then for the average propagator at time t one has the expression: 

oo n 

k(t) = ^2xn(t)Y[k v (At i ,V i ), 

n =0 i=0 


( 26 ) 
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where Xn(t) is the probability of having exactly n jolts up to time t. The above can 
be simplified if we recall our assumption that the time intervals and velocity values 
for each event are independent and identically distributed and the joint velocity-time 
distribution is [22] : 

P{v,t) = ^5(\v\ -uo)p(t), 

where p(t) is our heavy-tailed inter-event distribution function. Then one has: 


Kit) = Y J Xn{t)k v , 

n=0 



( K(At, v 0 ) + KjAt, - v 0 ) 

\ 2 


(27) 


where the averaging is now performed over the random time intervals, At, only. To 
proceed further we use the standard method in the theory of CTRW, namely work 
with time Laplace transformed quantities j22j. In the Laplace domain for example the 
probability Xn(t) has an especially simple form: y n (s) = p n (s)( 1 —p(s))/s [22]. This 
allows for writing the average propagator in the Laplace domain in a compact form: 


k(s) = jr l -jAL p n {s) X' = L_h£) (1 - p (s)kv)-\ (28) 

z — J s s 

n= 0 

where we have assumed that the geometric series of operators converges. Continuing 
borrowing from the established methods of the CTRW theory 1221 we note that the 
large time asymptote of the propagator is determined by small s behaviour of k(s). In 
particular by virtue of the so called Tauberian theorems the small s expansion of the 
fat tailed PDF p[t) oc T a /t 1+a is given by p(s) = 1 — r a s a + .... Then keeping only the 
principal terms in s one arrives at K[s) ~ T a s a ~ 1 ( 1 — K v )~ l which back in the time 
domain corresponds to the asymptote 

W)= (^“pW-jfi-i:)- 1 , H®. (29) 


Finally let us simplify the expression for the average single jump kernel K v in the moving 
basis (1251) . From the definition of K v in Eq. (127]) one obtains: 

where &k,k";v 0 is the characteristic function of the waiting time distribution, i.e. p(s), 
evaluated at the imaginary argument s* = {i/h){E{k) + E(k ) — E{k")). Assuming 
that we are interested only in the vicinity of the main six peaks (see above), that 
k ~ ko and additionally rE(ko)/h 1 one can use the familiar small s-expansion of 
the characteristic function and obtain: 

(k\l—K^\k') = -r a J dk"s°:{k,k"-v 0 ) [{k\^ x \k''){k''\e- iltx \k') + (k\e~ fkx \k"){k"\e i ^ x \k')^ (30) 

Note that according to Eq.(l29|) the average propagator decays in time as t~ a . It is 
interesting to compare it to the results of Section [4] for the Levy-Smirnov distribution 
with a = 1/2 as the average survival amplitude considered there is proportional to the 
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diagonal matrix element (ko\K (t)\k 0 ). One can see both from Eq. (l2T]) and Fig]4]that 
the survival amplitude indeed decays in time as t _1//2 as predicted by Eq. (|29lh What 
the latter formula fails to predict though are the oscillations of the corresponding Bessel 
function. To get these oscillations one must keep higher terms in the s-expansion of 
(1 — p(s)K v )~ 1 in Eq. (l28lh A similar situation in fact occurs when considering the large 
time limit of the characteristic function of a Levy-Smirnov walk (see Chapter 8 in [22]). 


5.2. The Brownian motion 


To get the ensemble averaged propagator K(k,k']t) (in the co-moving representation) 
for the case of Brownian motion ([2]) we first derive an exact evolutionary equation 
for the average amplitude c(k,t ) without assuming anything about the shape of the 
potential. This is done by averaging Eq.Q directly and applying Furutsu-Novikov 
formula (corresponding to Stratonovich regularization) [37]. The result reads (see 


Appendix C): 


i(k,t) = D / dk'(k\A 2 (t=0)\k')e i < Efk) - E(k ' w '‘c(k',t). 


(31) 


For completeness we have also derived an equation for the binary correlation function 
C(k, k']t) = c*{k',t) c(k,t) (see Appendix C ): 

— = d (A 2 (t ) c + cA 2 (t ) - 2 A(t)cA(t)) = d \A(t), \A(t), 6] ] , (32) 

at \ / L ' L ' JJ 

where we have used skew-Hermicity of A{t). This equation is quite interesting in itself 
but its detailed analysis lies beyond the scope of this paper. 

Let us use ( 1311 ) to calculate the the desired average propagator K(k, k 0 ; t ) which is 
the solution of the above equation subject to the initial condition c(k, 0 ) = 5{k — k 0 ). 
In the case of a delta-potential the matrix elements of the operator A 2 (t = 0 ) = d 2 jdx 2 
in the basis | k) can be computed analytically 

(k\A 2 (t = 0 ) 1 / 0 ') = — k 2 5{k — k') + (fl/n)t* k \ t\ k t\. ( 33 ) 

The fastest way to obtain this is to use the stationary Schrodinger equation for ijjk'{x) 
to get its second derivative and then project it onto ^ k ( x )• 

Then from Eq. flHTTl one immediately obtains an exact equation for the average 
coefficients in the moving frame: 


c(k, t ) = — D k 2 c(k, t) + D 


n 


IT 


dk 


, \kk , \e i ( E ( k )- E ( k ')'* t t n 

- tn)(\k'\ + iil) 


c(k', t). 


( 34 ) 


Given the initial condition c(k, 0) = S(k — ko) the solution can be presented as 


K(k , fc 0 ; t) = c(k , t) = S(k - k 0 ) e~ Dk2 A + 


DQ\k\ 


7T 


— itt) 


e -Dk\t-t’) e iE(k)t’/h g ^ dt ^ 
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where g{t) is the solution of the Volterra equation of the second kind: 


9(t) = 


f W\ e 


/1 p—i E(k')t/h 


U e ~iE(ko)t/h 


\k’\ 




Q — Dkgt 


dk' 


DQk' 2 

7T (k 12 + n 2 ) 


df'e- Dk ' 2 ( t - t ')- i ( E ( k 'y h )( t - t '') g(t'y 


(36) 


The hrst term in (1351) dehnes the average survival amplitude, f(k Q ,t ), and one 
immediately sees that it agrees with the BP result (ETT1) and the numerics FigJU The 
second term in (1351) is the averaged non-adiabatic correction (due to the diffusion in 
the fc-space) and in order to define its asymptotic behaviour one needs to know the 
behaviour of g{t). 

The Laplace time-transform is the method of choice and after some simple algebra 
one gets for the transformed function g(s): 

ko ko(l — i)DQ 2 


d(s) = 


(k 0 + + iO(k 0 )) + + iu(k 0 )) [v/2sa;(f]) + (1 + i)E(Q)/h 


,(37) 


where we have introduced the complex frequency oj(k) = E{k)/h — i Dk 2 . The time- 
domain function g(t) can then be obtained by the standard Mellin formula for the inverse 
Laplace transform. There is a simple pole at s = —iuj(k 0 ) in the left half plane. The 
contribution from this pole is the oscillating time decaying exponential: 

k 0 n (£ - w) 


0 i (<) = e 




(38) 


+(&~ iO) feofi]' 

A second contribution is obtained from the branch cut that is situated in the complex 
plane at s = —a/ (2u;(f2)) with a 6 R + . The sum of the two integrals over the two sides 
of the brunch cut provides the remaining contribution: 


02 (*) = 


fc 0 (l + i)DQ 2 


dxe 


x 


x~ 


V r(k Q + ill) tu(n)u;(/co) + 

It is possible to evaluate this integral in terms of the error function 
-fc 0 (l + i)Dtl 2 


x 2 - 2iuj(n)uj(k 0 ) x 2 + 2 


02(f) = 


g(Sl) 


2(/cq + iQ) (cu(h2)tu(/co) H —^2 


V 7 —2iu)(r2)d;(/co)e _ii: ’ (fco ^Erfc ( y/-itu{k 


+ V2ltA>e‘ ES l n >'/<'‘ 2 “( n >>Erfc ( 

n \ y n z u(\l) I 


(40) 


If we plug both gi(t) and (/2(f) into the expression (j35j) for the average propagator we 
obtain correspondingly the contribution ci(k,t) from g\(t) 

IjM e ~Dk$t e (i/h)(E(k)-E(ko))t 


Ci(k, t ) cx 


\k\ — ifl uj(k) — &(ko) 


.(39) 
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and the contribution c 2 (k,t) from g 2 (t) 


c 2 (k, t) oc 


\k\ 

\k\ — itt 


e iE 2 (p.)t/(h 2 u(Cl))+iE{k)t/h 

h 



+ v / ^(fij^)e-“ {fco)i+lE( ' :)t/s 


Erfc —iCj(ko)t^j 
(oj(k) - Cj(k 0 )) 


(41) 


The functions Ci(k,t) and c 2 (k,t ) are the corrections to the term localized at ko in Eq. 
(135|) . Both decay with time, ci(k,t) exponentially and c 2 (k,t ) oc 1/Vt for large t. 

To summarize, Eqs. (1291) and (135|) provide the analytical expressions for the disorder- 
averaged quantum propagator (in the co-moving representation) for the Levy and BM 
case, correspondingly. These expressions can be used e.g. to study quantum scattering 
of narrow wavepackets at a randomly moving scatterer which is much less studied area 
than the deterministic wavepacket scattering [131117]. 


6. The fluctuations of the pumped charge 

Let us now Anally look at the fluctuations of the pumped charge (J7]). According to 
the results of Section [4] the BP contribution of each scattering state \k 0 Cl leads 
to a static probability flux. On the other hand due to time-reversal symmetry of the 
static quantum dot potential the absence of bias makes the flow of the electric current 
impossible without applying an external gate voltage. Therefore it is the non-adiabatic 
correction to the amplitudes given by Eq. (1T81) that breaks the time reversal symmetry 
and introduces a random bias at each realization of the CTRW which in turn leads to 
the fluctuations of the pumped current and charge. For the unbiased CTRW the average 
values of both current and charge in ([7]) should vanish and this is indeed confirmed by 
the results of the numerical simulations (not shown). Therefore we shall be interested 
in the time growth of the fluctuations of the pumped charge, namely its RMS value, <tq. 
This value was sampled at boundary points ±L/2 chosen to be at least three times the 
standard deviation of the CTRW to avoid crossing effects. Since we have not noticed 
significant variation of the results for right and left lead (which is to be expected because 
of the current conservation) we only provide the results for the right lead x = +L/2 
which are given in Fig. [TJ 

We see from Fig {T]( a) that the time dependence of the standard deviation of the 
pumped charge shows a few common features for both the BM and Levy case. In 
particular, one observes that after a short transient time up to t ~ T F charge fluctuations 
seem to enter a stable growth phase with the slope larger than predicted by the snow¬ 
plow formula ([6]) i.e. 1/2 for BM and and close to it for the ballistic Levy case. In the 
latter case however the ballistic spread velocity is still almost an order of magnitude 
larger than predicted by Eq.([6]). We have been able to trace this regime for almost 
two decades of time variation. Of course the observed deviation from the BTP formula 
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Figure 7. The r.m.s of the pumped charge vs time measured in Fermi wave periods 
Tp = 2irh/Ep averaged over 4000 realizations, (a) Non-adiabatic regime, (b) adiabatic 
regime. 


should come as no surprise since the parameters of the simulation used in FigJT^a) were 
the same as in the previous chapters and were far from the adiabatic regime. To check 
what happens in the adiabatic limit we have repeated the simulation for a very slow 
moving quantum dot there the r.m.s. instantaneous wave vector k was much smaller 
than kp. The results are given in FigJTfb). There for the Levy case the effective wave 
vector was chosen to be k = 10 ~ 3 kp while for the Brownian motion we defined the 
instantaneous velocity via a typical r.m.s. displacement v = \J2D/t where r is the 
effective correlation length imposed by a finite time grid. In our dimensionless units the 
results of Fig0(b) are given for r = 4 x 1CT 3 , D = 10 -4 . As expected, asymptotically 
the snow-plow effect of Eq. (j6]) is recovered. Note that the transient period is roughly 
the same as in the non-adiabatic case of FigE^a). In both non-adiabatic and adiabatic 
cases this transient period corresponds to the time needed for the pump switched at 
t — 0 to start working and reflects the inertia of the pumping. In fact in the adiabatic 
regime this transient behaviour was already numerically observed in Ref.flT] for slow 
harmonic oscillations of the barrier strength and position. 

As for the large time asymptotes, we observed in the last chapter that the 
average propagator and the probability current seem to be determined solely by the 
BP contribution. Despite that, the charge fluctuations (that by virtue of expressions 
(p3|) , (f7T) and flU) are determined by the fourth order correlation functions of the reduced 
amplitudes c±(k,t ) = c(k, t) exp(± ikj(t))) are non-negligible and yield the charge 
fluctuations that grow faster than expected from the snow-plow formula (16|) . The 
exact theory of the phenomenon will require developing a kinetic description of the 
sums (integrals) of the fourth order correlation functions built on the c±(k,t) which will 
present some theoretical challenges due to the non-Gaussianity of their statistics. This 
however is beyond the scope of the present paper. 
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7. Conclusions and outlook 

In this paper we have studied the tunneling of a quantum particle through a randomly 
moving quantum dot and have shown that in the large time limit the disorder-averaged 
quantum state of a system behaves largely as if one simply applies the corresponding 
random geometric phase as is the case for the adiabatic motion of the scatterer. We 
have also studied the growth of the variance of the the pumped quantum charge (171) 
and found significant deviations from the “snow-plow” contribution ([6]) derived from an 
adiabatic BTP formula. We note that most of the results of this paper should apply 
to scatterers of an arbitrary (but well localized) shape provided that the typical spatial 
extent of the latter, A, is such that the inequality kpA -C 1 holds. While we only report 
here the results for the two random models Brownian motion and Levy walk we note 
that the results of the Chapter |4] (the separation of the BP contribution and the fast 
non-adiabatic fluctuations) are valid for an arbitrary random motion of the scatterer. 
Therefore we still expect that the average survival amplitude is asymptotically given by 
the characteristic function of the random walk, although the specific shape of it is of 
course process-dependent (see Eq. (l2lT) and FigH]). 

Concerning possible experimental realizations and especially the observation of the 
nonadiabatic Berry phase it is interesting to consider the analogy with optics, more 
precisely the propagation of an optical beam in a narrow randomly curved waveguide 
[38] where an equation of the form (JT]) is valid with time variable t replaced by the 
paraxial propagation distance z. Another interesting possibility arises from the unitary 
mapping of the considered problem of a randomly moving quantum dot to that of 
scattering at a spatially fixed barrier exposed to an external temporally fluctuating 
potential via a Kramers-Henneberger type transformation Hi. As for our results for 
quantum pumping of electrical charge, note that we get for the BM and Levy model for 
the values of parameters chosen rather high values of the charge fluctuations up to a 
few dozens of e. We think that such fluctuations should be observable experimentally 
in quantum dot experiments. 

Our work can be extended in various directions on the theoretical side: It would be 
highly desirable to quantify better the corrections determined by Eq. (fT8|) to the Berry 
phase result for the wavefunction and by doing so understand better the time dependence 
of the pumped charge fluctuations exceeding the predictions of the adiabatic theory. 
Furthermore testing our results for extended scattering potentials and considering the 
pumping of spin instead of charge would be interesting. 

We are grateful to Uzy Smilansky for initially drawing our attention to the problem 
of random tunneling and to Nick Korabel for useful comments on the Levy walk model. 
SD was supported by the Marie Curie Fellowship (project “INDIGO”) and wishes to 
acknowledge kind hospitality of the Faculty of Physics of the University of Duisburg- 
Essen. DW acknowledges the Minerva foundation for a fellowship during the time the 
project was initiated. 
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Appendix A. The saddle point contribution to the wavefunction 


The result in Eq. (fTfil) is calculated from c(k,t) by performing the /e-integration along 
the contour shown in Fig. [2] under the assumption that corrections collected along the 
path away from the real axis are zero. We want to determine here the contribution of 
these corrections for large times t. In this case we consider a saddle point approximation 
around k = 0. 

First one notices that the contributions from the two arches can be made 
exponentially vanishing by increasing the radius to infinity and one is left only with 
the contribution from the bisector: 


Sxkofat) = 


dk 


1 — i 


c(k(l - i)/V2,t)t k{1 _ iyV ^e 


ik(l—i){x— r y(t))/y/2 — hk 2 t/ (2m) 


(A.l) 


dk 


( 1-0 

V2 


c(ifc(l - i)/y/2,t ) (V fc(1 - i)(3; - 7(t))/v/2 + 7-J 1 _ i)fc/V 2 e “ ifc(1-i)( * _7(t))/ ' /2 ) 


with c(k,t) given by the last three lines in Eq. (1141) . We now expand the prefactor of 
e -hk 2 t/( 2 m) arounf ] ^} ie stationary point k — 0. The first two orders in k yield zero, up 
to second order we obtain for k > 0 

(1 _ j)c((l - i)k/V 2, « 

k 2 (1 + i) [Ar 0 (e** 07(f) - 1) - 2Qk 0 'y(t)d('y(t)) - 21} sin(koj(t))0(— 'y(i))] 

2 y/2 nk^fl^ko + ifi) 

and for k < 0 


(1 - i)c((l - i)k/V2, t) + r 

h 2 (1 + Q 2 [x — 7(0 + 1/(212)] . , . {t) 

2 V2 7rfc§(fcb + *n) ^ 

-2fisin(A;o7(O)0(-7(O)] • 


(l-i)jfc/\/2 fc 


- t —ik(l—i)(x—'y(t))/\/2 


/y/2 


- l) + 2Qk 0 'y(t)9(-'y(t)) 


Using the identity 


dkk\- hkH '^ = ^L(^\ 


3/2 


the final result for Sxk Q (x,t) reads 


4 y/2 \ht ) 7rkl(k() + iO) 


{ [fc 0 (e lkol{t) - l) - 2fl sin(/u O 7(t))0(—7(t))] 


[x ~ 7(0 + 1 /^] - 0)7(0 [0(7(0) “ 0(~7(O)] + 2Qkoj(t)0(-j(t))(x - 7 ( 0 )} 


(A.2) 


The expression above for Sxk 0 {x,t) is the leading order correction to the wavefunction 
in t resulting from the ^-integration for obtaining Xk 0 { x ,t) from c(k,t). Higher orders 
in 1/i are obtained by expanding c(k,t ) to higher orders in k. 

Concerning the probability flux (J8jl. we can evaluate the leading order current 
correction in t 


AJ ko = -A 
m 


d d 

X* ko (x, t ) Q^$Xk 0 {x, t ) + 5x* ko (■ x, t ) — Xk 0 [x, t ) 


dx' 


(A.3) 










25 


with Xk 0 {x, t) denoting the wavefunction obtained within Berry phase approximation 
in (JUD - These contributions are oscillatory in dependence of x and maximally increase 
with x like (xsin/co^)- 


Appendix B. The derivation of the “snow plow” formula for a slowly 
moving scatterer 

Let us show that for a barrier moving with constant small velocity the asymptotic 
result (|24|) reproduces correctly the adiabatic BTP formula (in particular the “snow¬ 
plow” result by Avron et. ah CD- According to Eq. (JTJ) one needs to integrate the 
probability flux of “left” and “right” scattering states under the ID Fermi surface to 
get the electric current at a given point. For a “left” scattering state with fc 0 > 0 we 
have already obtained an asymptotic formula (l24]h In order to find the “right” state 
we again use time reversal transformation i.e. the fact that if ij}(x,t‘,v) is a solution of 
the time-dependent Schrodinger equation (0) with an even potential then -0(— x, f; —v) 
is also a solution. At t — 0 and negative k 0 one has Xk Q {x, 0) = X\k 0 \(—x). Therefore 
the “right” scattering state can be obtained via its “left” counterpart (|24l) by replacing 
k 0 —> \k 0 \, x —> —X, k —> — k: 


Xk 0 {x,t ) = 


+ _ p -i\k 0 \x -iE(k 0 )t/h 

l \k 0 \+k e e 


x — y —oo 


(B.l) 


sjy | g i\ko\x e iE(ko,)t/h j ~^01^ iE(^ko\-\-‘2ki s )tll% ^ ^ j qq 

for t —> oo and k 0 < 0. Assuming that the pumped charge is sampled at a distant right 
lead x — L —> +oo and plugging (12^1 and (IB.ID into the equation for the charge (0) one 
gets: 

kp 

ik 0 \t ko _ k \ 2 + *( - ^o + (fc 0 + 2k)\r ko+k \ 2 + 


dQ 

dt 


eh 


2n m 


A 


2kr k j~Q 2i ( k O+k)L ^i(E(ka)-E(ko+2k))t/h ) 


dkn. 


When we consider the adiabatic limit where k <C fcp, one can expand the integrand 
leaving only the terms linear in k. Then one has (recall that \tk\ 2 + \rk\ 2 — 1 ): 

kp 


dQ 

dt 


eh 

2irm 


k 0 \t ko \ 2 -L, 0 ' k 0 k + /cq( 1 + |^fc 0 | 2 ) + 2/c|ry 0 | 2 


dk 


k + 2~m (r ko e 2ik ° L ) 


dko = 


ehk 


urn 


kp 


k F \r kr \ 2 + SR / 'r ht e 2 ' t ° L dk 0 


(B-2) 


The integral term in the last equation decays as 1/L and thus becomes negligible for L 
larger than a few multiples of the Fermi wavelength 2 Ti/kp- 

Finally, after an infinitesimally small time interval At the scatterer moves the 
distance Ay = ( hk/m)dt and from Eq. (IB.2D one gets exactly the snow-plow formula 



















Appendix C. The average amplitude and the correlation function for the 
Brownian motion 
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We start from the derivation of Eq. (l3Tf . From Eq.Q it follows that 


c(M) = €(t) 


J dk'A(k,k';t) c(k',t). 


(C.l) 


To average this equation over noise one can use the celebrated Furutsu-Novikov (FN) 
formula [37] which for the white Gaussian noise £(£) has the form: 

_ £ TD 

f(i)*ra = 2D—, (C.2) 

where /?[£] is an arbitrary functional of the noise. To apply this formula we need to 
calculate the functional derivative of the solution with respect to noise: 

s -j&± = linty_rt = lin!*/_>* 9[t -t'] [/ dk'A(k, kif) c{k £')] 

= | f A(k,k';t ) c(k',t)dk', 

where we have used causality and defined the value of the Heaviside function 0[O] = 1/2 
which corresponds to Stratonovich convention. Applying the above to the r.h.s. of (1C. II) 
Eq. ffdTTl follows. 

To get an equation for the correlation function we can follow a similar procedure 
as above. First we multiply the l.h.s of Eq. (|9l) by c*(k',t) add the hermitian conjugate 
and average over noise. We obtain: 


C(k,k']t) — / dk"A(k,k";t)$,(t)c(k",t)<f(k',t) + h.c.. 


Using FN formula (1C.2|) and the variational derivatives Sc(k,t)/5£(t) obtained above 
one arrives at Eq.( 
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